Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 24 February 2013 (MN IM^ style file v2.2) 



(N 



(N 



A Combined Optical and X-ray Study of Unobscured Type 
1 AGN. II. Relation Between X-ray Emission and Optical 
Spectra 

Chichuan Jin*, Martin Ward, Chris Done 

Department of Physics, University of Durham, South Road, Durham, DHl 3LE, UK 
Submitted to MNRAS 



o 



> 

m 
(N 

p 

rn 
O 
(N 



X 



ABSTRACT 

In this second paper in a series of three, we study the properties of the various emis- 
sion features and underlying continuum in the optical spectra of Type 1 active galactic 
nuclei (AGNs) by using the unobscured hard X-ray emission as a diagnostic. We in- 
troduce the use of the 'Correlation Spectrum Technique' (CST) for the first time. We 
use this to show the strength of correlation between the hard X-ray luminosity and 
each wavelength of the optical spectrum. This shows that for Broad Line Seyfert Is 
all the strong emission lines (broad component of Ha and H/3, [Nelll] AA3869/3967, 
[01] AA6300/6364, [Oil] AA3726/3729, [OIII] AA4959/5007) and the optical underly- 
ing continuum all strongly correlate with the hard X-ray emission. By contrast, the 
Narrow Line Seyfert Is show stronger correlation of the optical continuum but weaker 
correlation in the lines. 

A cross-correlation with luminosity between the various Balmer line components 
and the broadband SED components shows that the best correlation exists between 
the hard X-ray component and broad component (BC) of the Balmer lines. Such a cor- 
relation is weaker for the intermediate (IC) and narrow components, which supports 
the view that broad line region (BLR) has the closest link with the AGN's compact 
X-ray emission. The equivalent widths of Balmer line IC and BC are found to corre- 
late with L2-iokeV: ^2-iokeV ~ ^boi / L2-iokeV , Balmer line FWHM and black hole 
mass. There is a non-linear dependence of the Balmer line IC and BC luminosities 
with L2_iofcey and L5100, which suggests that a second-order factor such as the ILR 
and BLR covering factors affect the Balmer line component luminosities. The Balmer 
decrement is found to decrease from ~5 in the line core to ~2 in the extended wings, 
with mean decrements of 2.1 in BLR and 4.8 in ILR. This suggests different physical 
conditions in these regions, such as variations in electron density and dust abundance. 

The [OIII] line is composed of a narrow core together with a blue-shifted compo- 
nent with average outflow velocity of 130;tgo° ^^^ ^~^- "^^^ total luminosity of [OIII] 
A5007 shows the best correlation with the luminosity of hard X-ray emission, and so 
can be used to estimate the intrinsic X-ray luminosity of obscured AGNs. We use the 
CST to show the correlation of the [OIII] A5007 luminosity with each wavelength of 
the full continuum SED. This shows as before that the [OIII] A5007 strongly corre- 
lates with power law tail, but also that it correlates almost as strongly with the optical 
continuum from the disc, but not with the soft excess. 

Key words: accretion, Optical/X-ray correlations, active-galaxies: nuclei 



1 INTRODUCTION 

Active Galactic Nuclei (AGNs) are characterized by strong 
energy output over a very broad frequency range, from the 
infrared to hard X-ray. Although their spectral energy dis- 
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tribution (SED) has been studied for several decades, one 
of the main difficulties is the unobservable region between 
0.01 keV and 0.2 keV, which is due to Galactic and intrinsic 
photoelectric absorption. As a result most previous studies 
have focused mainly either on the UV /optical/infrared side 
or the soft/hard X-ray side, of the unobservable region. 
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1.1 Previous Work 

In the optical band, information about the AGN can be de- 
rived from the Balmer emission line series e.g. Ha, H/3 and 
H7, strong forbidden lines e.g. the [OIII]A4959/5007 dou- 
blets, and the optical underlying continuum. Timing analy- 
sis of variations in the Balmer lines via Reverberation Map- 
ping study is the most reliable method to estimate AGN 
black hole masses (e.g. Kaspi et al. 2000; Peterson et al. 
2004). Measurement of the [OIII] A5007 emission line is used 
to provide an orientation independent estimate of the AGN's 
central ionizing flux (Kauffmann et al. 2003; Heckman et al. 
2004; Brinchmann et al. 2004; Heckman et al. 2005). It is 
also a principal diagnostic of the properties of the narrow 
line region (NLR), such as the stellar velocity dispersion 
(e.g. Nelson & Whittle 1995; Nelson & Whittle 1996; Boro- 
son 2003; Greene & Ho 2005) and outflow speed (e.g. Bian, 
Yuan & Zhao 2005; Komossa et al. 2008). The underlying 
continuum in the optical is believed to be dominated by 
emission from the accretion disc (Koratkar & Blaes 1999; 
Jin et al. 2011, hereafter: Paper-I) and AGN's host galaxy 
(e.g. Hao et al. 2010; Landt et al. 2011). 

The soft to hard X-ray emission is emitted from a com- 
pact region. The hard X-ray power law tail above 2 keV is 
produced by disc photons which are Compton up-scattered 
by a hot, optically thin electron population (e.g. Haardt & 
Maraschi 1991; Zdziarski, Poutanen & Johnson 2000). How- 
ever, the origin of the soft X-ray excess below 1 keV is still 
not clear (e.g. Gierlihski & Done 2004; Done et al. 2011; also 
see Paper-I). Significant absorption of the continuum by ion- 
ized gas located near to the central region (i.e. the 'Warm 
Absorber') can further complicate interpretation of the ob- 
served X-ray spectra (Reynolds 1997; Crenshaw & Kraemer 
1999). 

The infrared band especially beyond 2 /im is believed 
to be dominated by emission from the dusty torus, and so 
provides information about its physical condition (e.g. Rees 
et al. 1969; Rieke 1978; Lebofsky & Reike 1980). The ob- 
served 1 lira minimum in the SED is also a common feature 
(Sanders et al. 1989). Recently, Landt et al. (2011) showed 
that the continuum luminosity at 1 ^m can be used to esti- 
mate the AGN black hole mass. 

Many previous studies have attempted to construct 
AGN broadband SED by combining multi-waveband obser- 
vational data to model the emission across the UV-soft X-ray 
energy gap (see references in Paper-I). Whichever method 
is used to construct the broadband SED, it is found that 
in almost all cases the SED peaks in the unobservable re- 
gion (rest-frame) which encompasses the 'Big Blue Bump 
(BBB)'. One important consequence of this is that the un- 
observable region contains the largest fraction of an AGN's 
radiated energy (Walter &: Fink 1993; Grupe et al. 1998, 
1999; Paper-I). 



1.2 The Approach Adopted in This Paper 

In Paper-I, we presented modelling results for a sample of 51 
unobscured nearby Type 1 AGNs. From the spectral fitting 
we derived numerous spectral parameters that were not con- 
taminated by X-ray or optical obscuration. In this paper, 
we investigate the link between the central ioniza- 
tion flux characterized by the unobscured hard X- 



ray, and the properties of the various optical emis- 
sion lines and the continuum. We conduct this inves- 
tigation by studying the profiles of various emission 
lines, and the luminosity correlations betAveen the 
hard X-ray continuum and various optical emission 
features. 

Narrow Line Seyfert Is (NLSls) are often considered as 
a special type of AGN whose permitted line width is compa- 
rable to other forbidden lines, and their [OIII] A5007/H/3 flux 
ratio is lower than is typical of normal Seyfert Is (Shuder 
& Osterbrock 1981; Osterbrock & Pogge 1985). The con- 
ventional definition for a NLSl is a Seyfert 1 with H/3 
FWHM<2000 km s^^ and [OIII] A5007/H/? <3 (Goodrich 
1989). Previous studies focused on such kind of Seyfert 
galaxies and showed many intriguing characteristics. For 
example, NLSls were often found to have low black hole 
masses (e.g. Wang & Lu 2001; Bian & Zhao 2004; Zhou et 
al. 2006; Komossa & Xu 2007; Zhu, Zhang & Tang 2009, 
hereafter: Zhu09). These may be systematically lower 
than predicted by the M-a, relation which holds w^ell 
for BLSls (Grupe & Mathur 2004; Mathur & Grupe 
2005). The NLSls also have high Eddington ratios (Boro- 
son 2002; Komossa 2008; Paper-I). It is thus proposed that 
NLSl's central black hole may still be growing (e.g. Mathur, 
Kuraszkiewicz, & Czerny 2001; Komossa & Mathur 2001; 
Komossa 2008). In addition, NLSls have softer 2-10 keV 
spectra, lower 2-10 keV luminosities, higher aox values and 
more energetic BBB (see Paper-I and references therein). In 
this paper, we continue to pay special attention to the NLSl 
subset in our sample and show how they behave differently 
from other sources in the cross-correlation study. 

This paper is organized as follows. We first review some 
most important characteristics of the sample in Section 2, in 
order to emphasize that our study conducted in the follow- 
ing sections should be related to the most intrinsic proper- 
ties of AGN's bare core. Section 3 will present the 'Optical 
to X-ray Correlation Spectrum (OXCS)' based on our new 
'Correlation Spectrum Technique (GST)', from which vari- 
ous correlation features related to the hard X-ray luminos- 
ity are found in the optical spectrum. Section 4 will study 
the cross-correlation between different Balmer line compo- 
nents and broadband SED components. Section 5 will focus 
on correlations related to the Balmer line equivalent width 
(EW). Section 6 will study the physical properties of differ- 
ent Balmer emission line regions. [OIII] A5007 line's property 
and its correlation with different SED components are put in 
Section 7. Summary and conclusion will be made in Section 
8. Following Paper-I, flat universe model is adopted with the 
Hubble constant Ho = 72 km s"^ Mpc"\ VLm = 0.27 and 
^A = 0.73. 



2 THE SAMPLE AND THE SPECTRAL 
MODELLING 

2.1 Sample Selection 

The sample used in this paper is a nearby unobscured Type 
1 AGN sample derived from the cross-correlation of 2XMMi 
& SDSS DR7 catalogs. The main selection steps are 
listed below for completeness. A full source list and 
more detailed explanation of the sample selection 
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Figure 1. Examples of spectral fitting in Paper-I. Left panel shows the broadband SED fitting of PG 1115+407 which consists of a 
modified accretion disc (green dashed line), a soft X-ray Comptonisation (orange dotted line) for the soft X-ray excess and a hard X-ray 
Comptonisation (blue dash-dotted line) for the hard X-ray power law tail. Right panel shows the emission line fittings of RBS 1423 
around Ho and H/3. Blue solid lines represent different line components. 



can be found in Paper-I. 

(1) We searched the 2XMMi and SDSS DR7 cata- 
logs and identified 3342 extragalactic sources having 
both X-ray and optical spectra. 

(2) Within these sources, we selected those with Hfi 
in emission and redshift z < 0.4, so that both the Ha 
and H/3 emission lines are covered by the SDSS spec- 
tra. This assists Avith modelling of the Balmer lines 
(see Paper-I). This selection resulted in 802 unique 
X-ray sources. 

(3) Within this sample set, we identified 96 Type 1 
AGNs all with a minimum of 2000 counts in at least 
one of the three XMM-Newton EPIC cameras, to 
ensure high X-ray spectral quality. 

(4) We then excluded 23 sources whose H/3 line was 
modified due to high reddening, Ioav S/N or a data 
gap in the SDSS spectra. The remaining sample con- 
tains 73 AGNs. 

(5) For each of the 73 sources, a power law model 
w^as fitted to the rest-frame 2-10 keV X-ray spec- 
tra. The 16 objects Avith photon index uncertainties 
greater than 0.5 were thereby excluded, leaving 57 
Type 1 AGNs w^ith relatively w^ell constrained 2-10 
keV spectra. 

(6) A further 6 objects w^ere excluded because of the 
obvious signature of a warm absorber at ~0.7 keV. 
The final sample contains 51 AGNs, with 12 AGNs 
being classified as NLSl using the conventional def- 
inition (Goodrich 1989), while the others are all 
BLSls. Most objects in this sample are radio quiet, ex- 
cept for 3 sources that were reported as radio loud, i.e. 
PG 1004-1-130, RBS 0875 and PG 1512-f370. High quality 
XMM-Newton EPIC X-ray spectra and SDSS optical spec- 
tra are available for every source in this sample. In addition. 



simultaneous optical/UV photometric data from the XMM- 
Newton OM monitor arc available for 37 sources. 

We exclude PG 1004-|-130 from all correlations 
as it is a BAL quasar, so its X-ray fiux is likely to 
be heavily obscured even though it does not show 
clear evidence for absorption edges (Miller et al. 
2006). We also exclude Mrk 110 from the optical 
correlations (Sections 3, 4 and 5) as this source 
show^s strong optical variability (Kollatschny et al. 
2001; Kollatschny 2003) and the SDSS spectrum is 
very different from the (non-simultaneous) XMM- 
Newton OM data (see Paper-I). However, the [OIII] 
line luminosity does not change with the optical con- 
tinuum, so we include this object in the [OIII] versus 
broad band SED correlations in Section 7. 



2.2 Selection Bias 



Our sample distributes evenly within 0.031 < z < 
0.377 with a mean redshift < z >= 0.137^0:^73. Its 
selection is mainly based on both high quality opti- 
cal and X-ray spectra, so any AGN that was not 
detected, or only marginal detected, by SDSS or 
XMM-Newton w^ould not be included in the sample, 
i.e. sources with low mass accretion rate or strong 
obscuration were excluded. Further selection criteria are 
more related to specific spectral characteristics, such as ex- 
cluding objects with heavy optical reddening. Type 2 objects 
and X-ray warm absorber objects. Therefore it will be hard 
to directly estimate the bias due to these selection effects. 
However, we could compare our sample's general properties 
with previous samples. We found that our sample had very 
similar redshift, 2-10 keV luminosity and bolometric lumi- 
nosity distributions as Vasudevan & Fabian (2007) 's sample 
(hereafter: VF07), except that VF07's sample also includes 
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some extremely nearby and low X-ray luminosity sources 
such as NGC 4395, NGC 3227 and NGC 6814, which did 
not fall into our sample. The redshift distribution of 
our sample is similar to that of Grupe et al. (2010), 
except that their sample has a higher fraction of 
lower redshift sources (< z >= 0.112 ± 0.077). The se- 
lection effect regarding the broadband SED shape should be 
weak, which is because the broadband SEDs of our sample 
have shown a very strong shape diversity, with the intrinsic 
optical to X-ray spectral index aox ranges 1~2. It is true 
that any objects with extraordinary optical to X-ray lumi- 
nosity ratios would not be included in our sample. But such 
odd broadband SEDs are more likely due to optical or X-ray 
obscuration, rather than being intrinsic to AGN, thus they 
are no longer the type of unobscured sources we need in this 
sample. 



2.3 Major Sample Properties 

The most important characteristic of this particular sample 
is the high quality of both their optical and X-ray spectra. 
In the optical, none of these sources suffer strong dust red- 
dening, thus all sources have very clear optical underlying 
continuum superposed by a series of clear broad and narrow 
emission lines. In the X-ray our selection criteria have ex- 
cluded sources whose spectra from XMM-Newton have low 
signal to noise or contain strong warm absorber features (i.e. 
the absorption edge at ~0.7 keV, from combined absorption 
features from partially ionized Oxygen and Iron, see e.g. Lee 
et al. 2001; Turner et al. 2004). The rest of the sources all 
have high quality X-ray spectra which represent the emission 
from the AGN's core emission. 



2.4 The Spectral Modelling 

Another important characteristic is the availability of all 
important spectral parameters from optical to X-ray for the 
whole sample, which results from our thorough modelling 
of the multi-component Balmer lines, optical spectrum and 
broadband SED. 

In the emission line fitting (e.g. Figure [l] right panel), 
two Gaussian components were used to fit the [OIII] A5007 
line, i.e. a central component and a blue component. Then 
the whole profile of [OIII] A5G07, i.e. including both 
central and blue components, was used in the fitting 
of the narrow component of Balmer lines. Two addi- 
tional Gaussian profiles were included in fitting each of the 
Ha and H/3 lines, so that each line contains a narrow com- 
ponent (NC), an intermediate component (IC) and a broad 
component (BC). All other strong nearby emission lines, e.g. 
[Nil] AA6585/6548 doublets, Li A6708, [SII] AA6717/6734, 
are included by adding more Gaussian profiles into the whole 
model. Various constrains were set for these Gaussian com- 
ponents which are all described in Paper-I. 

In the broadband SED fitting (e.g. Figure [1] left panel), 
we made use of a new SED model {optxagn model in Xspec 
vl2: Done et al. 2011), which modifies the accretion disc 
emission (the green dashed line) by assuming a corona ra- 
dius within which all accretion disc emission is transferred 
to a soft X-ray Comptonisation component (the orange dot- 
ted line) to account for the observed soft X-ray excess, plus 



a hard X-ray Comptonisation component (the blue dash- 
dotted line) to model the hard X-ray power law tail. We 
rebuilt the broadband SED from the optical to hard X-ray 
by extrapolating the best-fit model over the unobservable 
UV/soft X-ray region, and then derived all the SED param- 
eters. Detailed descriptions of these spectral fitting can be 
found in Paper-I. 



3 THE OPTICAL TO X-RAY CORRELATION 
SPECTRUM (OXCS) 

3.1 The Motivation of OXCS 

The hard X-ray emission from AGN rises from Compton 
up-scattering of disc photons by a high temperature (100s 
of keV) electron population which forms a corona region lo- 
cated above the accretion disc. Although this hard X-ray 
emission may only contribute a small fraction of the total 
central ionizing flux (depending on the corona radius, see 
Paper-I), it is capable of penetrating deeply into the most 
dense gas regions near AGN's core, which emits some specific 
line species. Provided that the intrinsic hard X-ray emission 
is not heavily absorbed and the optical spectrum is not heav- 
ily reddened, we can identify those optical emission features 
closely linked with the high energy core emission by test- 
ing their correlation with the hard X-ray. As noted above, 
our sample was selected based on the unobscured nature 
of both the optical and X-ray spectra, so the hard X-ray 
luminosity can be used as a reliable diagnostic to investi- 
gate its relation to the optical spectral properties. Here we 
propose and test a new type of spectrum, the 'Optical to 
X-ray Correlation Spectrum (OXCS)' (see Figure 01. This 
is a direct extension of previous monochromatic lu- 
minosity correlation studies between the X-ray and 
optical e.g. 'L2kev vs. L = (e.g. Green et al. 2009; 
Lusso et al. 2010). 



3.2 Construction of OXCS 

The principle behind the OXCS is to cross-correlate the hard 
X-ray luminosity (here we choose the luminosity of 2-10 
keV: L2-iofcey ) with the monochromatic luminosities at each 
wavelength of the optical spectrum for the whole sample of 
objects. Then we plot the correlation coefficient against the 
wavelength, to see how the correlation changes with wave- 
length. For each source in the sample, we corrected the SDSS 
spectra for Galactic reddening and de-redshifted them to 
their rest frame. We define a standard optical spectral re- 
gion that is covered by the SDSS spectrum of every source 
(around 3700-6700A), and calculate the monochromatic lu- 
minosity at 1000 wavelengths distributed evenly across this 
spectral range. The Spearman's rank test was used to cross- 
correlate these monochromatic luminosities with L2-iofeeV, 
and so the Spearman's p coefficient was derived at each of 
the 1000 wavelengths. Figure [2] plots the Spearman's p coef- 
ficient against the wavelength for the 12 NLSl (red line) and 
37 BLSl (blue line). The wavelengths of some of the most 
prominent optical emission and absorption features in a typ- 
ical AGN spectrum are indicated in the plot. Note that the 
spectral coverage is not exactly the same for different OXCS 
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Figure 2. The OXCSs for our sample, the method of constructing them is described in Section 13.21 The red line is 
the OXCS for the 12 NLSls in our sample. The blue line is the OXCS for the 37 BLSls in our sample. Purple and 
cyan dotted lines indicate the wavelengths of some most prominent optical emission lines for a typical AGN, writh 
cyan lines indicating the weaker line of any doublets. Green dotted lines indicates the wavelengths of Mg b and Na D 
stellar absorption features. '(NC)' is the narrow component, w^hile '(BC)' is the broad component. The dashed region 
of Ha(BC) means that this region is not covered by the BLSl OXCS. 



subsets, because these subsets have slightly different redshift 
ranges. 

3.3 Correlation Features 

Within the limited wavelength range of OXCS, the un- 
derlying continuum correlation would not be expected to 
change significantly, and so it forms a basic correlation 
level that is relatively flat in OXCS. Superposed on this 
basic correlation continuum there are various emission and 
absorption-line- like features, which shows that those lines 
have stronger or weaker correlation with L2_iofeeV than 
the underlying continuum. We identify some of the most 
noticeable characteristics in the OXCSs as below: 

(1) \j2-wkev emission correlates well with the entire 
optical underlying continuum. The underlying cor- 
relation is ~0.8 in the BLSl OXCS, and ~0.9 

in the NLSl OXCS. A noticeable phenomena is that 
the underlying correlation does not decrease significantly 
towards either the blue or red end of the optical spec- 
tral range. This confirms that our sample suffers little 
intrinsic reddening, and that the optical continua redward 
of 5000A are dominated by AGN emission, i.e. the host 
galaxy contamination is small in most cases. There is also 
suggestions about the existence of an extra component 



contributing optical emission at redward of 5000A, which is 
probably originated from the self-gravity dominated region 
of accretion disc (Vanden Berk et al. 2001; Collin & Hure 
2001; Puchnarewicz et al. 2001; Soria & Puchnarewicz 2002; 
Pierens, Hure & Kawagushi 2003; Collin & Kawaguchi 2004; 
Hao et al. 2010; but see Landt et al. 2011). But we find it 
difficult to investigate this component merely using SDSS 
spectra due to the difficulties in accurately subtracting the 
host galaxy emission. 

(2) In the BLSl OXCS, the broad wings of Ha and 
H/3 correlate better with L2_iofceV than the optical 
continuum, and so result in the apparent broad-wing-like 
features around 4860A and 6565A. However, the core 
region of Balmer lines has a much weaker correlation with 
L2-iofceV as shown by the two narrow correlation dips 
centred at 4862.68A and 6564.6lA in the OXCSs. This 
directly shows that the Balmer lines in BLSls consist of 
(at least) two components, from the NLR and BLR of 
different physical conditions. We will investigate this issue 
further in later sections. The [OIII] AA4959/5007 doublets 
in BLSls show a very strong correlation with L2-iofceV, iu 
spite that the [OIII] AA4959/5007 originates in the NLR 
far from AGN's core. We will investigate this in more 
detail in Section [7] The Balmer line profile in NLSls is 
not strongly broadened, so the NLSl OXCS does not 
exhibit similar correlation bumps around Balmer 
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lines as seen in the BLSl OXCS. In contrast to the 
BLSls, the [OIII] AA4959/5007 in NLSls also have much 
weaker correlation with L2_iofeeV than their local optical 
continuum. 

(3) The BLSl OXCS also exhibit emission-line- 
like features at the wavelengths of some other emission 
lines in a typical BLSl optical spectrum, such as [Nelll] 
AA3869/3967, [Ol] AA6300/6364 and [Oil] AA3726/3729. 
This suggests that these emission lines all correlate strongly 
with the hard X-rays. However, [OI] A6300 is a relatively 
weak line, and its prominence in the OXCS may support the 
existence of dense gas clouds near AGN's core inside which 
gas stays neutral or at low ionization. Only hard X-rays can 
penetrate into these clouds and produce such low ionization 
lines. The fact that [Oil] AA3726/3729 correlates quite well 
with \j2-wke.v suggests that reddening is indeed quite low 
for our sample, since otherwise the presence of dust would 
tend to diminish any correlation. It is apparent that in 
terms of the OXCS around emission lines, the NLSls are 
different from BLSls. This may be a result of geometrical 
effects, or that the line emitting regions in NLSls are not as 
closely associated with hard X-ray emission as in the BLSls. 

(4) On the contrary, it is seen that the stellar absorp- 
tion lines Mg b and Na D do not correlate well with 
L2-iofeei/, producing absorption-like features in the OXCSs. 
Neither does Fell emission in the ranges 4400-4800 A and 
5100-5600 A correlate with L2-iofceV, especially for BLSls. 



3.4 The Correlation Spectrum Technique (CST) 

More generally, the OXCS provides a new tool for spectral 
studies based on medium to large samples of objects. We will 
name this the 'Correlation Spectrum Technique (CST)'. As 
shown above, an example of the CST is the OXCS, which 
has proved to be useful for investigating hard X-ray related 
correlations in the optical spectrum for different AGN pop- 
ulations. 

Using hard X-ray as a diagnostic, we can also apply the 
CST to spectra of longer wavelengths such as near and far 
infrared, provided that an AGN sample with good spectral 
data is available. We can also use luminosities other than 
hard X-ray as the diagnostic in CST. For example, in Sec- 
tionOwe constructed the 'SED to [OIII] A5007 Correlation 
Spectrum (SOCS)', in which case the CST uses the luminos- 
ity of [OIII] A5007 to produce the correlation spectrum from 
optical to hard X-ray for different sample subsets. 



4 BALMER LINE LUMINOSITY 

For Type 1 AGN, each Balmer line profile consists of two 
to three distinct components, i.e. narrow component (NC) 
from the narrow line region (NLR) which extends a few hun- 
dred parsecs from the central black hole (e.g. Bennert et al. 
2006; Heckman et al. 2005), broad component (BC) from 
the broad line region (BLR) which is tens of light-days from 
the black hole (e.g. Kaspi et al. 2005; Bentz et al. 2006). 
Sometimes another intermediate component (IC) is present, 
with moderate linewidth which is probably originated from 
the intermediate line region (ILR) which may extend up to 



the region of the inner radius of the dusty torus (e.g. Zhu09) . 
Therefore the distance of these line emitting regions from the 
compact AGN goes as NLR, ILR, BLR from the farthest to 
the closest. Among these emission line regions we may also 
expect a density gradient and a correlation trend with the 
central ionizing flux. The correlation between the Balmer 
line luminosity and X-ray luminosity has been known long 
ago (e.g. Ward et al. 1988), but the hard X-ray correla- 
tions for different Balmer line components have never been 
studied. Now our sample with both high quality optical and 
X-ray spectra provides an opportunity for this study. 



4.1 Balmer Line Component Luminosity vs. 

L2-10fceV 

We have shown in the previous section that different com- 
ponents in the Ha and H/3 lines may have different corre- 
lation status with L2_iofceV- In this section we investigate 
the correlation of L2-iofceV with the individual components 
(Narrow, Intermediate and Broad) derived from our Balmer 
line decompositions for each source in Paper-I. The best 
correlation is for L2-iofceV vs. H/3 IC-I-BC (i.e. the combi- 
nation of NC and BC, equivalent to NC subtracted H/3) 
luminosity, shown in the upper panel in Figure O NLSl 
(red squares) and BLSl (blue circles) clearly lie on the same 
strong correlation. Spearman's rank test gives a correlation 
coefficient of ps — 0.9 and probability of random distribu- 
tion of d, = 8.1 X 10~^^ 

We also plot the uncorrected data from our two ex- 
cluded sources (PC 1004+130: purple star and Mrk 110: 
green circle) on the correlation. These strongly deviate from 
the best-fit line (orange solid line), the ilcr lines (orange 
dash lines) and the ±2cr region (light gray region). How- 
ever, PG 1004-1-130 is the only BAL quasar in our 
sample, whose X-ray emission after correcting for in- 
trinsic absorption, w^as reported as 0.73 dex weaker 
than normal PG RLQs normalized to similar opti- 
cal/UV luminosities (Miller et al. 2006). Mrk llO's 
optical continuum is highly variable (Kollatschny et 
al. 2001; Kollatschny 2003). The SDSS spectrum of 
Mrk 110 is ~1 order of magnitude less luminous than 
the FAST optical spectrum (Landt et al. 2011), which 
is just the required amount of correction we need to pull 
Mrk 110 back onto the best-fit correlation line. Therefore, 
we conclude that this tight correlation between H/3 IC-f BC 
luminosity and 2-10 keV luminosity is an intrinsic property. 
We propose that if an AGN is found to strongly deviate 
from this correlation, then it is likely that its X-ray or opti- 
cal emission is obscured. 

We calculated three types of regression lines as de- 
scribed in Isobe et al. (1990). Since Ha also gives similar 
results (see Table [T]), we give bisector regression equations 
for both H/3 and Hq: 
(i) L2-iofceV expressed by the Balmer line Luminosities: 

L03L2-10 = {0.83±0m)LogLHo.(NCsub) + (8.35±1.43) (1) 
Lo(?L2-io = i0.83±Qm)LogLHp(NCsub) + (8.56±1.52) (2) 
(ii) Balmer line Luminosities expressed by L2_iofceV : 

LogLHc,{NCsub) = (1.20±0.05)LoffL2-io - (9.50±2.09) (3) 
LogL„fnj,csub) = (1.18±0.04)Lo3L2-io - (9.17±2.04) (4) 
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Figure 3. The luminosity correlation between Balmer line components and 2-10 keV. The upper figure shows H/9 luminosity (NC 
subtracted) vs. L2_iofceV- The connected filled and empty purple stars indicate the position of PG 1004+130 before and after being 
corrected for the 0.73 dex (Miller et al. 2006). The connected filled and empty green circles indicate different optical positions of Mrk 
110 as calculated from the SDSS spectrum and the FAST spectrum (Landt et al. 2011). The solid orange line shows the linear regression 
line treating L2_iofeeV a-s the independent variable, with the two dashed orange lines indicating the itlo" region for new observations, 
and the shaded region showing the ±2cr region. The lower panels present the same type of correlations for different H/9 components, i.e. 
H/3 NC, IC, BC, IC+BC (or NC-sub) and the whole line. In each plot. Spearman's rank coefficients were calculated after excluding PG 
1004+130 and Mrk 110. The regression coefficients are listed in Table [B] 



Using the above equations we estimate that PG 1004+130's 
X-ray luminosity is weaker than normal Type 1 AGNs 
by 1.0±0.3 dex, which is slightly higher than Miller et 
al. (2006) 's estimation of 0.73 dex weaker than normal 
PG radio-loud quasars with similar optical/UV luminosities 
based on their X-ray spectral analysis. 

The lower panels in Figure [3] show the (weaker but still 
very significant) correlations for the different line compo- 
nents. It is clear that the BC has the best correlation with 
hard X-ray emission as found previously; IC+BC and the 
whole H/3 line also show good correlations; the NC related 
correlations are not as good as others. This confirms that 



the BLR has the closest link with the AGN's central X-ray 
continuum. 



We note that the well-known Malmquist bias (Gonzalez 
& Faber 1997, and references therein) will be partly respon- 
sible for these correlations, but it should affect all of the cor- 
relations equally. Hence the change in correlation strengths 
among the different Balmer line components should be real. 
We also examined these correlations using flux instead of lu- 
minosity, and a very similar trend of strength of correlation 
was found for NC, IC and BC. Note that the components in 
Hq show similar correlation status as in H/3 (see Table [!}. 
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Table 1. The luminosity correlations between Balmer line com- 
ponents and 2-10 keV. The Bisector regression coefficients in the 
equation: l^og{LBalmerLine)=Po'^og{L2-iokev)+io are listed, 
along with the Spearman's rank correlation coefficients (ps, ds) 
as defined in Table [B] 



Line Comp. 
Ha 


Bisector Re 

/3o 


gress Coef. 
?0 


Rank C 

Ps 


or. 

ds 


NC vs. La-iofceV 


1.00±0.06 


-2.01±3.28 


0.76 


-10 


IC vs. L2_10feeV 


1.18±0.05 


-9.00±2.26 


0.92 


-20 


BC vs. L2_10feeV 


1.23±0.05 


-11.51±2.19 


0.93 


-22 


IC-I-BC vs. La-iofcev 
(i.e. NC-sub) 


1.20±0.05 


-9.50±2.09 


0.94 


-22 


Whole vs. L2-iofceV 


1.15±0.04 


-7.66±1.95 


0.93 


-22 


H/3 


/3o 


& 


ps 


ds 


NC vs. L2_10fceV 


1.02±0.07 


-3.49±3.35 


0.76 


-10 


IC vs. L2_10fceV 


1.17±0.05 


-9.41±2.54 


0.86 


-15 


BC vs. L2-10fceV 


1.21±0.05 


-10.67±2.09 


0.93 


-21 


IC-HBC vs. L2_l0fceV 

(i.e. NC-sub) 


1.18±0.04 


-9.17±2.04 


0.92 


-20 


Whole vs. h2-u)kev 


1.15±0.04 


-7.91±1.97 


0.93 


-21 



4.2 Cross-correlation between Balmer Line 
Components and Broadband SED 
Components 

In Paper-I we decomposed the Balmer lines into broad, in- 
termediate and narrow components. We also decomposed 
the broadband SED into three components, namely the disk, 
soft X-ray Comptonisation and hard X-ray Comptonisation. 
Thus it is possible for us to correlate each line component 
with each SED component. The results are presented in Ta- 
ble [B] and Figure lAll Note that for each source we 
subtracted the Fell emission from the nearby re- 
gion of the H/3 line before conducting the profile de- 
composition, so the dispersion in these correlations 
is a result of Fell contamination. We find that among 



the three SED components, the hard X-ray Comptonisa- 
tion produces the best correlations, while the accretion disc 
emission and soft X-ray Comptonisation show weaker cor- 
relations. Among the different Balmer line components, the 
correlation strengthens from NC, IC to BC. Figure |4] shows 
the correlation status for H/3 IC-I-BC vs. broadband SED 
components. It is clear that the best correlation is found 
in the hard X-ray component though the correlation is also 
good for adding both soft and hard X-ray Comptonisation 
together as total corona luminosity. 

FigureHalso shows that for BLSls the H/3 IC-hBC cor- 
relate well with the accretion disc luminosity and bolomet- 
ric luminosity. However, NLSls are much more dispersed 
in these correlation plots, and thus dilute the correlation 
strength of the whole sample. We need a larger sample to 
confirm the different behaviors of NLSls and BLSls in these 
correlations. 



5 BALMER LINE EQUIVALENT WIDTH (EW) 
5.1 Balmer Line Component EW vs. 'L2-iokeV 



In the previous section we reported the strong correlation 
between H/3 luminosity and L2_iofeeV as hnp oc LjL^^pj^y"* 
for the whole H/3 line, and ^hp(ic+bc) oc L2;^Jq^°,^°* for the 
IC plus BC in H/3. The index 1.18±0.04 is bigger than unity 
which indicates that the luminosity of broader components 
of H/3 increase faster than linearly with L2-iofcv (also see 
TablelTJ. The situation is the same for Ha. In order to inves- 
tigate this issue further, we study the properties of Balmer 
line EW. We perform the cross-correlation between L2-iofceV 
and H/J component EW. Figure [5] shows the results. There 
is no correlation between L2-iofceV and H/3 NC EW, but the 
correlations between L2-iofceV and H/3 IC, BC EWs are sig- 
nificant as confirmed by Spearman's rank test (see Table IB)) . 
A bisector regression analysis shows that H/J IC-I-BC EW 

T 0.45±0.03 
^ ^2-lOfeV' ■ 

We also find clear anti-correlations be- 
tw^een the 2-10 keV bolometric correction (i.e. 
«:2-iofceV-=L6oi/L2-iQfeeV') and the H/3 EWs, as shown 

is the fraction 



in the second row of Figure [SJ k 

of 2-10 keV emission in the Bolometric luminosity, thus 



2-lOkeV 
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Figure 5. The correlations of H/3 component EW vs. L2_iofceV (first row), K2_iofeeV (second row). Different symbols represent the 
same type of sources as in Figure [3] Spearman's rank coefficients arc calculated for the whole sample. The orange dotted line indicates 
the bisector regression line. 



these correlations suggest that as the fraction of 2-10 
keV emission increases, the EWs of H/3 IC and BC also 
increase. This is not surprising since there is no correlation 
between H/3 EWs and the bolometric luminosity, as shown 
in Table |B] the Spearman's rank correlation coefficients 
(p^) are only -0.08, 0.13, 0.21, 0.22 and 0.22 for U^; vs. 
the EW of H/3 NC, IC, BC, IC-^BC and the whole line, 
respectively. The sources with large K2-iofceV are mostly 
NLSls which tend to have high mass accretion rates in 
terms of Eddington, 'Lboij^Edd (Vasudevan & Fabian 2007, 
2009; Paper-I; Jin et al. in prep., hereafter: Paper-Ill). 
However, when we directly cross-correlated Eddington ratio 
with the various H/3 EWs, no significant correlations were 
found (Figure |A2|) . Similar results were found when using 
Ha line instead of H/3. 

5.2 Does A Balmer Line Baldwin Effect Exist? 

The existence of the H/3 Baldwin effect is controversial. It 
was reported by Zhu09 that the H/3 IC EW anti-correlates 
with the monochromatic luminosity at 5100A (hereafter: 
L5100), with Pearson rank correlation coefficient being -0.48 
and at 99% level of smaller than 0, while the BC EW and 
the whole NC subtracted H/3 EW does not correlate with 
Lsioo- They suggest this should be due to the fiat geome- 
try of ILR and spherical geometry of BLR, but they did not 
mention whether such anti-correlation could also be found in 
their Ha IC or not. To compare with Zhu09's results, we use 
our Balmer line fitting results to perform similar correlation 
test. 

We first investigate the correlation between L2-iofceV 
and L5100 which can be seen directly from the OXCSs in 
Figure [2I The L2-iofceV vs. L5100 correlation is plotted in 
Figure [S] with Spearman's ps = 0.88 and ds = 1.1 x 10"^". 
The bisector regression lines are found to be: 
(i) L5100 expressed by L2-iofeeV-: 



LogiLsioo) = {0.92±0.05)Log{L2-iokev) + (3.76±2.24)(5) 

(u) L2_iofceV expressed by L5100: 

Log{L2-iokev) = (1.08±0.06)Lo(?(L5ioo) ~ (4.07±2.66)(6) 

Considering the strong correlations reported in previous 
paragraphs between H/3 IC, BC EWs and L2-iofeeV, the cor- 
relations between H/3 IC, BC EW and L5100 were expected. 
However, we do not confirm any strong positive or negative 
correlations between H/3 IC, BC EWs and L5100, though 
there is a large scatter in these cross-correlation plots, as 
shown in Figure IA2I and Table |B] There is a weak anti- 
correlation between NC EW and Lsioo- We therefore con- 
clude that no evidence of Baldwin effect is found in our study 
for Balmer IC and BC, but there may be such an effect for 
NC. 

In fact, there are some important uncertainties in 
Balmer line decomposition and L5100, which need to be con- 
sidered before studying the Baldwin Effect of Balmer line 
components: 

(1) The spectral quality of the Balmer line profile is crucial 
since currently decomposition of Balmer line totally depends 
on the line profile. This is more of a problem for H/3 since 
strong reddening can significantly reduce its S/N and dis- 
tort its profile, so the H/3 decomposition for such reddened 
sources will be unreliable. But this is not a problem for our 
low reddening sample. 

(2) Even with high quality Balmer line profiles, the line de- 
composition still has uncertainties. It is highly probable that 
the Balmer lines in Type 1 AGN must contain at least a NC 
and a BC. Assuming a Gaussian or Lorentzian profile for 
the BC, an additional IC is required during the line fitting 
procedure by the x^ statistics. However, the assumption of 
Gaussian or Lorentzian profile is not secure. In some cases, 
Balmer lines also exhibit a double-peak profile (Eracleous & 
Halpern 2003; Strateva et al. 2006; Bian et al. 2007) or an 
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Figure 6. L2_iofceV vs Lsioo- Different symbols represent the 
same type of sources as in Figure |3] Solid orange line is the bi- 
sector regression line assuming L2_iofcel/ is the independent vari- 
able. In each histogram, the red region highlights the distribution 
of the 12 NLSls in our sample. 



extended flat red wing (e.g. Mrk 0926, see Paper-I), which 
causes problem for the three-component decomposition. 
(3) The L5100 may not have just one contributor. In addition 
to the standard accretion disc emission in the optical, L5100 
may also contain stellar emission from host galaxy. It was 
reported that an additional component, probably from the 
outer region of accretion disc where self-gravity dominates, 
might also contribute a significant fraction of L5100 (see Sec- 
tion E3J. 

Therefore, it is difficult to find the intrinsic correlations be- 
tween Balmer line component EWs and L5100 due to the 
above uncertainties, and so the existence of Baldwin effect 
in any Balmer line component is still unclear. 



5.3 Balmer Line Component EW vs. FWHM and 
BH Mass 

Another result we found for Balmer line EW is its cor- 
relation with Balmer line width and black hole mass. 
We cross-correlate the H/3 component EWs with the H/3 
FWHM/c+sc and the 'best-fit' black hole mass (see the de- 
scription of 'best-fit' black hole mass in Paper-I). Figure [7] 
shows our results. It is clear that there are significant cor- 
relations between FWHM, black hole mass and the EWs of 
IC and BC. The best correlation is again found in BC. The 
results suggest that as the black hole mass increases (so does 
the Balmer line width), the emission from ILR and BLR be- 
comes more luminous relative to the continuum luminosity. 
We also note from Figure[7]that if only consider BLSls (cir- 
cular points in Figure [7]), then these is almost no correlation 
either between H/3 FWHM and component EWs or between 
black hole mass and H/? component EWs. The broadest 12 
BLSls (blue circular symbols) exhibit larger scatter than the 
rest of the sources. However, the correlation between H/3 BC 



EW and black hole mass is very strong for NLSls (red square 
symbols). There also seems to be a weak anti-correlation be- 
tween black hole mass and H/3 NC EW. Similar results can 
be found if replacing H/3 with Ha. 



5.4 The Nonlinear Dependence of Balmer IC and 
BC Luminosities on L2-iofceV and L5100 

As shown in previous sections that the relations between 
Lif^(ic+BC), L2_iofceV and L5100 can be expressed as: 



JH"/9(/C+SC) 



rl.l8±0.04 rl-28±0.05 



(7) 



Similar results can be found for Hq. Such non-linear de- 
pendences imply that if the continuum luminosity is the 
first-order factor, then there must be a second-order fac- 
tor causing the EW of Balmer line IC and BC to depend 
on L2-iofeeV- This second-order factor could be the covering 
factor of the ILR and BLR seen from the central ionizing 
continuum. AGNs with higher L2_iofeeV and L5100 may also 
have larger ILR and BLR covering factors, making their 
Balmer IC and BC EWs larger than in the AGNs of low 
L2-iofceV and Lsioo- 



6 THE PROPERTIES OF ILR AND BLR 

6.1 Balmer Decrement 

In addition to the luminosity and EW of Balmer line IC 
and BC, we can also investigate the intrinsic properties of 
ILR and BLR by studying the profiles of Ha and H/3. The 
Balmer decrement is one of the main parameters 
to investigate. It was reported that the change in 
Balmer decrement may arise from changes in the 
physical conditions of the partially ionized line emit- 
ting regions (Kwan &: Krolik 1979; Kwan &: Kro- 
lik 1981; Mathews, Blumenthal & Grandi 1980 and 
Canfield & Puetter 1981). For example, a decrease 
of Balmer decrement may be due to an increase of 
electron density Ne~, or an increase of ionization 
parameter H (i.e. the ratio of the photon density 
to the gas density), or an increase of Lya optical 
depth TLya (Krolik & McKee 1978; Davidson & Net- 
zer 1979; Shuder82). A high Balmer decrement can 
also be explained by a high dust abundance. This ar- 
gument was used as an evidence in Zhu et al. (2009) 
to support the link betw^een ILR and dusty torus. 
Therefore, Balmer decrements can be used as a clue 
to infer the physical conditions of the emission line 
region, and so it is important to obtain an accurate 
measurement of the decrement value. We calculated 
the Balmer decrement between Ha and H/3 for each line 
component. Figure [8] shows our results. The five histograms 
from top to bottom show Balmer decrement (Ha/H/3) dis- 
tributions of NC, IC, BC, IC-HBC and the whole line. The 
mean Balmer decrements with 1 standard deviation are also 
listed in Table H 

The results show that although the Balmer decrements 
of the whole Balmer line distribute around 3, the situation 
for different line components is quite different. The NLR has 
a big range of Balmer decrement values, with most probable 
value lying between 3.5-4.5, and a mean value of 4.85±1.81. 
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Figure 7. The correlations of H/3 component EW vs. H/3 IC+BC FWHM (first row) and best-fit black hole mass (second row). Different 
symbols represent the same type of sources as in Figure[3] Spearman's rank coefficients were calculated for the whole sample. The orange 
dotted line indicates the bisector regression line. 



Table 2. The mean Balmer Decrements of different line com- 
ponents for NLSl population, BLSl population and the whole 
sample. 



Line Comp. 


NLSl 


BLSl 


Whole Sample 


NC 


4.31±1.09 


5.02±1.96 


4.85±1.81 


IC 


3.50±0.61 


5.24±2.33 


4.83±2.18 


BC 


2.19±0.52 


2.11±0.91 


2.13±0.84 


IC+BC 

(i.e. NC-sub) 


2.73±0.30 


3.12±0.79 


3.03±0.73 


Whole Line 


2.94±0.30 


3.16±0.72 


3.11±0.65 



Therefore, it may imply the presence of some dust in the 
NLR. The IC from ILR also has a big mean decrement 
of 4.83±2.18, while The BC from BLR has a small mean 
decrement of 2.13±0.84. However, the Balmer decrement of 
the IC+BC distributes around 3, which seems that the low 
decrement in BC and high decrement in IC are artificially 
due to our multi-Gaussian decomposition. To investigate 
this issue, we conduct the following study. 

First, we think that the Balmer decrement we found for 
NLR is correct, based on the fact that this component has 
a narrow width, and is matched to the observed [QUI] pro- 
file. For BLSls, the profile of the NC is easy to define as a 
small spike superposed on the much broader line wings (see 
the Balmer line profiles in Paper-I); while for NLSls, untan- 
gling of the NC may introduce bigger uncertainties. How- 
ever, since the NLR is far from AGN's centre and thus much 
less sensitive to the AGN's central properties, we should not 
expect big difference in NLR's properties between BLSls 
and NLSls, such as Balmer decrement. This is consistent 
with what we see in Figure [S] that the mean decrement val- 



ues are similar between NLSl and BLSl. This supports our 
assertion that the NC decomposition is reliable, and the de- 
rived high Balmer decrement in NLR is real for the majority 
of sample sources. 

Second, we have reason to conclude that our decrement 
distributions for the IC and BC are also intrinsic to Balmer 
line regions. To prove this, we must first note that our com- 
bined H/3 and Ha line fitting has ensured very similar line 
decompositions for the two Balmer lines. To be specific, the 
IC and BC have the same central velocity shifts and FWHM 
in both H/3 and Ha, but the IC and BC can have very differ- 
ent relative fluxes (for more detailed description see Paper- 
I) . The observed decrement distribution differences between 
IC and BC reveal real changes in the decrement values across 
the emission line proflle. 

To see this directly, we divided each of the two Balmer 
lines into 10 segments in the velocity space from -5000 km 
s~^ to 5000 km s~^, and calculated the decrement value 
in each segment. We performed this spectral analysis for 
each object in our sample. Then the average decrement value 
in each segment was calculated for sources with a reliable 
decrement measurement in that segment. This method is 
model-independent except for the subtraction of the local 
underlying continuum and removal of the [Nil] A6584,6550 
doublets which uses the line fitting results from Paper-I. The 
mean values for NLSl subset and BLSl subset in each seg- 
ment are plotted in Figure |9] We see that the decrement 
peaks at the line centre and then decreases towards both 
sides. This suggests a low decrement in the broad wings, 
which is mainly modelled by the BC, and so supports the 
low decrement value found for the BC. For BLSl subset, 
the decrement in the red side (indicating a positive veloc- 
ity) is lower than in the blue side (negative velocity). It 
also appears in Figure [^ that NLSls tend to have lower 
Balmer decrements than BLSls in each segment. However, 
Figure |8] shows that the BC decrement distributions for the 
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Figure 9. Balmer decrement changing across the Balmer line 
profile from + 5000 km s^^ to — 5000 km s~^ . Each data point 
represents the average decrement value in that segment with the 
vertical bar showring the ±1 standard error. Blue points show the 
results for the BLSls; red points show the results for NLSls. But 
due to the small line width of NLSl, the Hux outside +/— 3000 
km s^^ for NLSls is of low S/N, thus only the mean decrement 
values in the central +/— 3000 km s~^ region were calculated 
and shown. The horizontal purple dotted line is a reference line 
at F(a)/F(/3) = 3. 



Figure 8. Balmer decrement distributions of different Balmer 
line components. In each panel the entire histogram shows the 
distribution of the whole sample, with green solid line indicating 
the mean decrement value. Red region highlights the distribution 
of the 12 NLSls and the orange dashed line indicates their mean 
decrement value. The cyan dashed line shows the mean decre- 
ment value of the BLSls. The mean decrements are also listed in 
Tabled 



NLSls and BLSls are similar, although they are sometimes 
of low contrast to the continuum. Thus it suggests that the 
observed lower Balmer decrement of NLSls is mainly due 
to the lower decrement value of the IC in NLSls, which is 
also consistent with the IC decrement distributions in Fig- 
ure[5] Therefore, we conclude that the observed differences in 
Balmer decrement distributions between different line com- 
ponents are mainly due to the complex decrement status 
across the Balmer line profile, with the broad wings having 
a lower decrement value and contributing mainly to the BC. 
There is also support for our findings from previous 
work by Shuder (1982) (hereafter: Shuder82), who showed 
ratios of the Ha and H^ lines for 18 Seyfert 1 AGNs. They 
reported that the average Hq/H/3 ratio ranged from 4.8 in 



the core to 2.2 in the wings, which is very similar to what we 
have found. They explained this as approaching the inner re- 
gion of AGN, the velocity dispersion would increase, and so 
do the Ne~ and H. This is also consistent with the system- 
atic inflow velocity we found for the BC as we will discuss 
later. Zhu09 also reported that the Balmer decrement in the 
IC was 4.78 and in BC it was 2.54, although their line de- 
composition for Ha and H/3 were not linked as in our study. 
However, they explained this high decrement in IC as due 
to the higher dust reddening in the ILR, and so could sup- 
port the link between ILR and dusty torus. Therefore, the 
Balmer decrement change can be explained by either dust 
abundance or line optical depth processes. Since it is likely 
that BLR is closer to the core region than the ILR, timing 
analysis such as detailed reverberation mapping of Balmer 
lines can be used to distinguish the radii between ILR and 
BLR. If the ILR is connected with BLR, then the changes 
of physical parameters such Ne~ , H and TLya rnay be the 
explanation of the higher Balmer decrement in ILR than 
in BLR. Alternatively, if the ILR is confirmed to be a dis- 
tinct region from BLR and is close to the dusty torus, then 
dust reddening may be a more plausible answer. Finally, we 
find no correlation between the Balmer decrement and the 
bolometric luminosity, as was suspected by Shuder82. 
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Table 3. The fraction of each Balmer line component and its 
correlations with black hole mass (TsAbh)} bolometric luminosity 
(Ljoi) and the H/3 IC+BC FWHM. ps and ds are Spearman's 
rank coefficients as explained in Table iBl 





frac 


vs. M 


BH 


vs. 


^bol 


vs. FWHM 


Hq 




ps 


ds 


Ps 


ds 


Ps 


ds 


NC 


10±8 


-0.6 


-6 


-0.2 


-1 


-0.7 


-7 


IC 


53+11 


0.1 


-0 


-0.1 


-0 


0.1 


-0 


BC 


37+11 


0.4 


-2 


0.3 


_2 


0.4 


-2 


IC+BC 


90+8 


0.6 


-6 


0.2 


-1 


0.7 


-7 


H/3 


% 


ps 


ds 


ps 


ds 


Ps 


ds 


NC 


7+6 


-0.6 


-5 


-0.2 


-1 


-0.5 


-4 


IC 


38+12 


-0.0 


-0 


0.0 


-0 


-0.4 


-2 


BC 


55+14 


0.4 


-2 


0.1 


-0 


0.5 


-5 


IC+BC 


93+6 


0.6 


-5 


0.2 


-1 


0.5 


-4 



6.2 Balmer Line Component Fraction 

Zhu09 proposed an evolutionary scenario for the emission 
line region, which claimed that as the black hole mass and 
luminosity increased, ILR and BLR would become closer 
to each other and finally merge. In our sample 49 of the 
51 sources show the necessity of including two broad Gaus- 
sian components to fit their Balmer lines, as suggested by 
the Bayesian Information Criteria (BIC). The fraction of 
each Gaussian component is calculated for both Ha and H/3, 
and cross-correlated against the black hole mass, bolomet- 
ric luminosity and the component line width. The results 
are listed in Table[3l The dominant component in H/3 is IC, 
but in Ha it is BC. The NC fraction is ~10% in both lines. 
We confirm a weak anti-correlation between the IC 
fraction and IC FWHM of H/3 (Spearman's p, = -0.4, 
ds — —2), similar to that found by Zhu09. But no such 
anti-correlation was found in Ha. So we conclude that it is 
still not clear whether the anti-correlation between the IC 
fraction and IC FWHM is an intrinsic property of the ILR. 
We did not find any correlation in the IC or BC fractions 
vs. black hole mass or bolometric luminosity. Therefore, the 
scenario proposed by Zhu09 regarding the geometry of ILR 
and BLR can not be confirmed by the results of our study. 
Instead, our results suggest that the ILR may simply be 
an intermediate region between BLR and NLR regardless of 
the black hole mass and bolometric luminosity. As the black 
hole mass increases, the luminosities of NC, IC and BC in 
both Ha and H/3 all increase, but the luminosity of IC and 
BC increase more significantly than the NC, resulting in the 
significant anti-correlation between the black hole mass and 
NC fraction. 



H/3 lines with more extended wings, along with higher Ed- 
dington ratios; while sources in the second population have 
broader but fiat-topped H/3 lines, together with lower Ed- 
dington ratios. They found a weak anti-correlation between 
FWHM/cTKne and Eddington ratio where aune is the second 
moment of H/3 line (see the definition given in Peterson et 
al. 2004). This can be explained in terms of higher turbu- 
lence in the core region of high Eddington ratio AGNs. We 
can explore this result for our sample. 

The second moment was measured from the NC 
subtracted H/3 line profile. It is clear from the defini- 
tion that the flux farther from the line centre would 
have more contribution to the total second moment, thus 
o";ine/FWHM/c+BC can be a representative of the broad 
wing strength compared to the whole line profile, i.e. a 
higher value of aune/FWHMic+BC corresponds to stronger 
broad wings. In our Balmer line profile fitting, the broad 
wing is mainly modelled by the BC, thus the FWHM 
ratio between BC and IC should have similar physi- 
cal meaning as cr;i„e/FWHM/c+BC- We correlate both 
aH„JFWRMic+BC and FWHMisc/FWHMic with Ed- 
dington ratio (Lioi/l^Edd)- 

Figure [10] shows our results. Orange points are the 
binned values over X-axis (Eddington ratio) with one stan- 
dard error in the Y-axis (line width ratio). The errors in 
Eddington ratio and line width ratio for each data 
point are dominated by systematical uncertainties 
in the spectral fitting which is difficult to quan- 
tify (see Paper-I), therefore individual error-bars 
are not shown. Such uncertainties will increase the 
dispersion in this correlation plot. However, positive 
correlations have been confirmed in both panels with Spear- 
man's rank test being ps =0.47, ds = 10~^ (upper panel) 
and Pi,=0.35, ds = 10~^ (lower panel). Our results confirm 
the correlation between H/? line profile and hboi/^Edd- In- 
terestingly, we also find similar correlations when replacing 
Fboi/^Edd with K2-iokeV and aox- For aii„e/FWiiMic+BC , 
the Spearman's test has ps=0.46 {ds = 10~^) for K2-wkeV 
and ps=0.51 [ds = 10"*) for a^^. For FWHMsc/FWHM/c, 
the Spearman's test has ps=0.54 {ds = 10""*) for K2-iokeV 
and ps=0.54 [ds — lO"**) for aox- But since the Eddington 
ratio is an intrinsic AGN parameter, and it also correlates 
with both K2-iokeV and aox (Vasudevan & Fabian 2007; 
Vasudevan &: Fabian 2009; Lusso et al. 2010; Grupe 
et al. 2010), it may be the driving parameter that regulates 
the Balmer line shape, as was also suggested by Collin06. 
Therefore, the FWHM of the BC may depend primarily on 
black hole mass, but it may also be regulated by \jboi/^Edd- 
This also explains the stronger correlation between black 
hole mass and the IC FWHM than between black hole mass 
and the BC FWHM. 



6.3 Balmer Line Shape 

6.3.1 Dependences on Eddington ratio 

The complex Balmer line profile can often be well modelled 
using contributions from the NLR, ILR and BLR, but it 
is also likely that local turbulence in the BLR may further 
broaden the Balmer line profile and be responsible for the 
presence of a very broad wing. Collin et al. (2006) (hereafter: 
CoUinOB) divided the reverberation-mapped sample into two 
populations: sources in the first population have narrower 



6.3.2 Inflow Implied by the Balmer Line Profile 

Another effect that may change the width of the whole 
Balmer line is the systematic velocity structure of the BLR 
as evinced by the general redshift of the BC and IC. It was 
found previously that both the IC and BC may be associ- 
ated with inflows (e.g. Sulentic et al. 2000; Hu et al. 2008). 
Our analysis also shows that for the whole sample both IC 
and BC have a wide range of velocity shifts relative to the 
central component of [OIII] A5007, but on average we find a 
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statistically significant shifts of 100 km s~^ for the IC and 
550 km &^^ for the BC. The inflow velocity we flnd for the 
BC is also consistent with the ~400 km s~^ typical inflow 
velocity of the Fell emission features found by Hu et al. 
(2008), supporting their conclusion that Fell emitting re- 
gion may trace some portions of the BLR exhibiting inflow. 
It seems probable that there is a velocity gradient within 
the BLR clouds, with the inner region of the BLR having 
a higher inflow speed, which gives rise to the extended red 
wing. However, the multi- Gaussian Balmer line decomposi- 
tion method used in this study cannot resolve the detailed 
changes in the kinematics and physical conditions within 
BLR. Therefore, a much more detailed broad line spectral 
and timing study is required. The relative velocity shifts 
between the IC and BC determines the asymmetry of the 
Balmer lines, but we did not find any significant correla- 
tion between the relative velocity shift and Eddington ratio. 
So the relative velocity shift cannot explain the correlation 
between aune/FWHMic+BC and hboi/'^Edd- 



6.4 Summary of results for the ILR and BLR 

To summarize the previous subsections, we propose the fol- 
lowing characteristics of the ILR and BLR. First, these two 
regions are likely to be closely related, with the ILR being an 
extension of the BLR. We found no evidence to support the 
ILR to be a distinct region from BLR. The inner region of 
the BLR may produce the red wing of the Balmer profile, in- 
dicating an systematic infiow velocity. The physical param- 
eters change continuously from the ILR into the BLR, prob- 
ably with increasing electron density, ionization parameter 
and Lya optical depth. The infiow velocity of the ILR gas is 
smaller than the BLR. The possibility that the ILR is associ- 
ated with the dusty torus cannot be ruled out, and indeed its 
Balmer Decrement is higher than found for the BLR. Unfor- 
tunately we are unable to draw any firm conclusions about 
the geometry of the ILR or BLR. A second-order factor such 
as the covering factor of the ILR and BLR may cause the 
faster than linear dependence of the Balmer IC and BC lumi- 
nosities on the continuum luminosity. Considering the tight 
correlations between the IC FWHM, BC FWHM and black 
hole mass, both the ILR and BLR should be gravitational 
bound and Virialized. For the BC from BLR, its FWHM 
may further be affected by the Eddington ratio through the 
process of local turbulence. Different inflow velocities of the 
ILR and BLR may also modify the shape of the Balmer 
lines. 



7 PROPERTIES OF EMISSION LINE [OIII] 
A5007 

The NLR may extend hundreds of parsecs from the AGN's 
compact core. It is ionized by the central continuum in a 
bi-conical geometry with an axis defined by the plane of the 
dusty torus. Since the NLR extends far from the dusty torus, 
the intrinsic dust reddening is expected to be low. Therefore 
in both Type 1 and Type 2 AGNs, the luminosity of narrow 
optical emission lines from NLR can provide an orientation- 
independent estimate of the central ionizing radiation (e.g. 
Mulchaey et al. 1994; Heckman 1995). 
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Figure 10. The H/3 line shape correlation with Eddington ra- 
tio. The upper panel uses FWHMse/FWHM/o to represent H/3 
shape, while the lower panel uses oune/^^^^lC+BC instead. 
In each panel the various symbols represent the same type of 
sources as in Figure |3] The orange data points are the binned 
data for different Eddington ratio bins with 1 standard error on 
the Y-axis. 



7.1 [OIII] A5007 Component Luminosity vs. 

L2-10fceV 

As one of the strongest narrow forbidden lines, [OIII] A5007 
is often employed as a proxy to estimate the intrinsic lumi- 
nosity of type 2 AGN (e.g. Heckman et al. 2004; Brinchmann 
et al. 2004). This is not only because of its large EW, but also 
because it is free from serious contamination of other spec- 
tral features. Heckman et al. (2005) (hereafter: Heckman05) 
used [OIII] A5007 as an optical selection criteria to study the 
difference between optical and X-ray selected AGN samples. 
They showed a tight correlation between [OIII] A5007 lumi- 
nosity and hard X-ray luminosity for Type 1 AGNs in both 
optical and X-ray selected samples. But the correlation is 
much weaker for an optically selected Type 2 AGN sam- 
ple, which is mainly due to their X-ray weakness resulting 
from intrinsic photoelectric absorption. Since our Type 1 
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Figure 11. The luminosity correlations between components of [OIII] A5007 and h2~iokeV- AH symbols and lines have the same 
meanings as in Figure |3] In each panel, a histogram is shown for the Log(LrQjjn;^5Q07/L2_iofcel') values of our sample, with the red 
histogram highlighting the NLSls. 



AGN sample has been carefully selected based on the high 
quality optical and X-ray spectra and the absence of se- 
vere absorption, our [OIII] A5007 vs. rest frame 2-10 keV 
luminosity correlation should be indicative of intrinsic con- 
nections. We also find that the line profile of [OIII] A5007 
consists of two components, a dominant central component 
and a blue-shifted component. Thus we analyze these two 
components separately. 

Figure [11] shows our results. We find that both com- 
ponents in [OIII] A5007 have very strong correlations with 
^2-iokeV, as confirmed by the Spearman rank test. The to- 
tal [OIII] A5007 line correlates slightly better with L2-iofc<:i' 
than the two separate components, suggesting that the cen- 
tral ionizing radiation ionizes both the outflowing [OIII] 
A5007 region (which produces the blue component), and the 
spatially more extended [OIII] A5007 region (which produces 
the central component). The solid orange line is the regres- 
sion line assuming L2-iofte\^ to be the independent variable. 
The two dashed orange line indicate the ±lcr region for new 
observation. The shaded region denotes the ±2cr region. 

Again, we put our two excluded objects on the correla- 
tion. Mrk 110 (green circle) sits on the best fit line as nei- 
ther its X-ray nor its [OIII] line luminosity are affected by 
the optical continuum variability. The X-ray weakness of PG 
1004-1-130 (purple star) does not cause strong deviation in 
these correlation plots considering the dispersion. The other 
outlier is IRXS J122019 (the blue circle point farthest from 
the shaded region) whose SDSS spectrum shows that this is 
one of the 'broadest' BLSls, but its narrow lines (including 
[OIII] A5007,4959 doublets) are much weaker relative to the 
broad lines than any other sources in our sample. It is pos- 
sible that compared with other sources, IRXS J122Q19 has 
a smaller NLR covering factor. 

We also calculated Log(L[o///]A5007/L2-iofcev) using 
the luminosity of the whole [OIII] A5007 line. We de- 
rived a mean value of -1.88±0.31 for the whole sam- 
ple, and -1.78±0.30 for the 12 NLSls. This is consistent 
with but slightly lower than -1.59±0.48 reported by Heck- 



man05 based on their sample of 20 Seyfert Is, which is 
likely due to the fact that their sample includes sources 
with strong X-ray absorption, e.g. they included Type 1 
AGNs such as NGC 3227 and Mrk 766 whose X-ray spec- 
tra are absorbed by warm absorber. Trouille & Barger 
(2010) (hereafter: TBIO) reported -1.85±0.5 for their 19 
BLAGN sources with z < 0.5 and -1.76±0.5 for their 
~100 BLAGNs with z < 0.85. Georgantopoulos & Aky- 
las (2010) (hereafter: GAIO) reported -1.98±0.39 for their 
34 Seyfert Is. These values are all consistent with ours 
within la. Our mean Log(L[oi/i]A5007/L3-20fcev) value is 
-2.00±0.33, which is consistent with -1.96'-^2.14 reported by 
Heckman05. Note that compared with previous works, our 
sample is spectrally 'cleaner', i.e. more carefully selected 
for low reddening and absorption sources. Thus our mean 

Log(L[O///]A5007/L2~10fcev) and Log(L[o///]A5007/L3_20fcev) 

are better constrained as having smaller dispersions. 

The correlations between [OIII] A5007 and L2-iofceV can 
be used to estimate the intrinsic hard X-ray luminosity es- 
pecially for nearby Type 2 AGNs and calculate the X-ray 
luminosity function (e.g. Sazonov & Revnitsev 2004; Shi- 
nozaki et al. 2006; Yencho et al. 2009; GAIO). We present 
our OLS bisector regression lines below: 
(i) L[o/J/]A5007) expressed by L2-iokeV and I/3-20fcev: 

LogLioiii]5007 = (1.06±0.05)Lo3L2-io - (4.44±2.64) (8) 
LogLioni]5oo7 = {1.02±0.05)LogL3-20 - (2.93±2.65) (9) 
(n) L2-iokeV and -L3-20feeV expressed by L[o///]a5007): 
LogL2-w = {0.94±0. 05) Log Lioni]5007 + (4.20±2.26) (10) 
LogL3.20 = (0.98±0.05)Lo3i[o///i5007 + (2.87±2.43) (11) 

Note that since the luminosity measurements of both [OIII] 
A5007 and L2_iofceV niay contain uncertainties from the in- 
trinsic absorptions and variability, the OLS bisector regres- 
sion method is more appropriate than the standard OLS 
method used in previous works. Combining previous re- 
sults from Heckman05, TBIO and GAIO with ours, we can 
conclude that such tight luminosity correlations between 
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Figure 12. The 'SED to [OIII] A5007 Correlation Spectra (SOCS)'. This is produced by calculating the Spearman's rank coefRcient 
between the [OIII] A5007 luminosity and the luminosity contained in each energy bin of broadband SED, thus the bigger coefficient 
indicates the better correlation in that energy bin. Lines of different color show the SOCS of different subsets as been labelled in 
the plot. The r2_iofceV^2.0 subset (SI: red line) contains 16 AGNs; the r2_iofeeV*Sl-8 subset contains (S2: blue line) 18 AGNs; the 
1.8 < T'2-iokeV < 2.0 subset (S3: orange line) contains 16 AGNs. Only spectral ranges below 0.006 keV and above 0.3 keV have 
observational data. The ionizing flux responsible for [OIII] A5007 emission is above 0.035 keV as shown by the purple dotted line. The 
two shaded regions are where model extrapolation was used. 



hard X-ray and [OIII] appear valid for at least z < 0.85, 
Log(L[O///]A5007)=38~44 and Log(L2-iofeev)=40~46. 

7.2 The SED to [OIII] A5007 Correlation Spectra 

(SOCS) 

The ionizing energy for [OIII] A5007 is 0.035 keV, thus all 
photons above this energy can in principle produce [OIII] 
A5007 emission. As another application of CST, we cross- 
correlate [OIII] A5007 luminosity with the luminosity con- 
tained in the continuum SED in each energy bin from optical 
to X-ray, and produce the 'SED to [OIII] A5007 correlation 
spectrum (SOCS)'. Our broadband SED model consists of 
three components: the disc, soft X-ray Comptonisation and 
hard X-ray Comptonisation. The soft and hard Comptoni- 
sation components can both contain a significant amount of 
energy above 0.035 keV, the disc emission may also extend 
above this energy when the black hole mass is small and 
mass accretion rate is high (Done et al. 2011). We divided 
our sample into three subsets: r2-iofeev^2.0 (SI: 16 AGNs), 
r2-iofc.v<1.8 (S2: 18 AGNs) and 1.8 < Ta-iofc.v < 2.0 (S3: 



16 AGNs). PC 1004-1-130 is excluded due to its unique SED 
shape but Mrk 110 is now included since the previous section 
shows that its [OIII] and hard X-ray luminosity are not dis- 
torted by the optical continuum variability. SI includes all 
12 NLSls, while S3 contains the broadest BLSls. A SOCS 
was calculated for each of the three subsets. 

Figure [12] shows the resultant SOCSs. Note: in this 
study only spectral ranges below 0.006 keV and above 0.3 
keV have observational data. Overall, [OIII] A5007 correlates 
best with hard X-rays above 2 keV. It is also well-correlated 
with the optical emission. However, the correlation in the 
UV/X-ray region is poor, which may be caused by the spec- 
tral modification due to Galactic and intrinsic extinction. 
This indicates that the hard X-ray ionizing flux also has a 
strong link with the optical flux which is presumably domi- 
nated by accretion disc emission. 

Regarding the SOCS of different subsets, we find that 
SI has the strongest correlations in the optical/UV band, 
which implies that our broadband SED fitting may be more 
reliable for sources in SI whose soft X-ray excess is more 
likely to be a real extra component (Middleton et al. 2009; 
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Jin et al. 2009; Middleton, Uttley & Done). The S2 group 
shows highly significant correlations in the hard X-ray band, 
which may imply that the hard X-ray power law tail of 
r2-iofceV^l-8 is an intrinsic separate component rather than 
being an artifact caused by absorption or reflection (Done et 
al. 2011). Therefore, the SI and S2 groups may indeed rep- 
resent two distinct types of AGNs (e.g. NLSls and BLSls). 
We also find that the correlation for S3 in optical and X-ray 
bandpasses is much less significant than for either SI or S2. 
This may indicate that other spectral factors such as absorp- 
tion and reflection, may be more important for the sources 
in S3, in which case our three-component SED model is too 
simple to recover their intrinsic SEDs. 



7.3 Outflow of NLR Implied by [OIII] A5007 
Profile 

The two components in the profile of [OIII] A5007 have been 
reported previously. Bian, Yuan & Zhao (2005) indicated 
that these two components are related to two physically 
distinct regions. Komossa et al. (2008) also reported blue 
outliers in [OIII] A5007 whose blue-shift velocity is up to 
SOO'^-^IOOO km s~^, favoring a decelerating wind NLR sce- 
nario. These results are all confirmed by our study. We find 
that for our sample, the velocity shift of the blue com- 
ponent in [OIII] A5007 relative to the central component 
ranges from -610 to -0 km s~^, and the mean velocity is 
— 130^230 ^''^ ^~^ ■ W^ ^^'^ ^"^^ ^ strong correlation be- 
tween the FWHM and velocity shift of the blue component, 
as Spearman rank test gives: ps=0.52 and da = lO"''. The 
larger FWHM of the blue component implies a smaller dis- 
tance from AGN's core region, so this correlation suggests 
that the outfiow velocity of inner NLR emitting [OIII] A5007 
is higher than that in the outer NLR. An outfiow speed de- 
creasing as it fiows away from the centre is a signature of 
decelerating wind. 



8 SUMMARY AND CONCLUSIONS 

In this paper, we made use of the detailed spectral fitting of 
an AGN sample reported in Paper-I, to study their optical 
spectral properties using their hard X-ray luminosity as a 
diagnostic. Our study focused on the H/?, Hq and [OIII] 
A5007 emission lines and the underlying continuum. The 
main results are summarized below. 

• The OXCSs have been constructed for different sub- 
sets of AGNs using our new spectral analyzing technique 
called CST. The OXCSs reveal many correlation features 
with L2-iofceV across the entire optical spectrum. Some were 
known previously, others are new. For example, the en- 
tire optical underlying continuum strongly correlates with 
L2-iofcev. [Nelll] AA3869/3967, [01] AA6300/6364, [Oil] 
AA3726/3729, [OIII] AA4959/5007 and the IG and EG in 
Balmer lines all well correlates with \j2-iokeV especially 
for BLSls. However, stellar absorption lines. Fell and the 
NC in the Balmer lines has much weaker or no correlation 
with L2-iofceV. We find some evidences for differences in the 
OXCSs between NLSls and BLSls. 

• The cross-correlation between luminosities of H/3 and 
Ha line components and the broadband SED components 
were performed. The results suggest that among the three 



SED components, the hard X-ray power law component cor- 
relates the best with Balmer line luminosity, and the corre- 
lations strengthen from the NC, IG to BG of Balmer lines. 
This supports the view that the BC has the closest link with 
AGN's central UV /X-ray continuum emission. 

• Significant correlations were found between the H/? 
component EWs and L2-iokeV, 1^2-iokeV' H/? FWHM and 
black hole mass, although these correlations become weaker 
for the BLSl subset alone. By cross-correlating Balmer line 
component EWs with L5100, no evidence for the 'Baldwin 
Effect' was found for the IG and BG, but such effect is weakly 
detected for the NG. 

• Our results suggest a faster than linear dependence 
of Balmer line IG and BG luminosities on the underlying 
continuum (e.g. L2-iofceV and L5100; Equation[7]), implying 
the presence of a second-order factor. We propose that this 
second-order effect could be the covering factor of the BLR 
and ILR seen by the central UV/X-ray continuum, so that 
higher L2_iofceV and L5100 sources may also have larger ILR 
and BLR covering factors. 

• We carried out detailed Balmer line shape studies in 
order to reveal the nature of ILR and BLR. We found that 
the Balmer Decrement value, defined by Hq/H/?, peaks at 
the line centre and decreases towards both sides, with the red 
wing having a lower decrement than the blue wing for BLSl 
subset. This was also consistent with IG's average decre- 
ment value of 4.83±2.18 compared to the BC's 2.13±0.84. 
These results, along with the systematic infiow speed we 
found in the BC (mean velocity: 550 km s~^), support the 
scenario that the inner region of BLR forms the red wing 
while the outer edge links with the ILR. Compared to the 
ILR, the BLR may have higher infiow speed, higher electron 
density, larger ionization parameter or higher Lya optical 
depth, A weak correlation between the shape of Balmer line 
profile and Eddington ratio was confirmed. A higher Ed- 
dington ratio corresponds to a more extended wing relative 
to the overall Balmer line structure. This implies that the 
velocity width of the Balmer line is not simply determined 
by the black hole mass, but also affected by local turbulence 
whose strength depends on the Eddington ratio. The higher 
Balmer decrement in ILR than in BLR could also be ex- 
plained as ILR has higher dusty abundance, but we found 
no other evidence to support ILR's link with the dusty torus. 
A weak anti-correlation between the EW of Balmer line NG 
and black hole mass was found. 

• In our study of [OIII] A5007, we refined its tight corre- 
lation with h2-iQkeV and h3,-20kev- We found that the blue 
and central components of [OIII] A5007 should be added 
together to provide the best correlation with hard X-rays. 
Using our best-fit broadband SEDs from Paper-I, we pro- 
duced the SOGSs for different sample subsets. The SOGSs 
show strong correlations between [OIII] A5007 luminosity 
and the continuum luminosities in either optical or hard X- 
ray bandpass. Subset SI and S2 both have highly significant 
correlations in the hard X-ray band, which implies that the 
shape of hard X-ray power law tail in these two subsets are 
intrinsic in spite of their totally different photon indices. But 
the SED of moderate V2-wkeV sources in S3 may be more 
complex. 

• The mean outfiow velocity of the blue component in 
[OIII] A5007 is -130+^1^0 km s~\ The strong correlation be- 
tween the FWHM and velocity shift of the blue component 
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in [OIII] A5007 suggests that the outflow speed of [OIII] 
A5007 clouds decreases from the central region outwards, 
suggesting a decelerating wind. 

• In this paper, we present well constrained equations 
which can be used to convert between the luminosity of 
Balmer line broad component and the intrinsic h2-iok<iV 
(Equation 1~4), between the intrinsic L5100 and L2_iofceV 
(Equation 5~6), and between the [OIII] A5007 luminosity 
and the intrinsic L2-iofceV and \jz-20kev (Equation 8~11). 
We suggest that these equations be used for inferring the 
intrinsic optical and X-ray luminosities of obscured sources 
such as BAL quasars or Type 2 AGNs, and for calculating 
the X-ray luminosity function. Considering the limited red- 
shift range of and size of our sample, similar studies should 
be carried out on larger samples to test the robustness and 
evolution of these equations at high redshift, which requires 
high quality infrared spectra. 
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Figure Al. The luminosity correlations between H/3 line components and SED components. Red points represent NLSls; blue points 
represent the broadest H/3 line BLSls; green point is Mrk 110; purple star is PG 1004+130; purple symbols indicate radio loud sources. 
The orange dotted line denotes the OLS regression line assuming the SED component luminosity is the independent variable. Spearman's 
rank correction coefRcient ps for the whole sample is also given in each panel. 

APPENDIX A: SUPPLEMENT OF BALMER COMPONENT CORRELATION PLOTS 
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Figure A2. The cross-correlation plots between H/3 line component EWs and L5100A (the monochromatic luminosity at 5100A) and 
^/^Edd (the Eddington ratio) Different symbols have the same meaning as in Figure \KT\ Spearman's p is given in each panel. 
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APPENDIX B: THE SPEARMAN'S RANK CORRELATION MATRIX BETWEEN Ha, H/3, [OIII] A5007 

LINE COMPONENTS AND SED COMPONENTS 
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Table Bl. Spearman's Rank Correlation Coefficients of SED components vs. Ha, H/3 And [OIII]5007 Line Parameters, ps is Spearman's rank correlation coefficient with 1 being perfect correlation and -1 
being perfect anti-correlation, ds is the logarithm of the significance level of being random distribution, i.e. a smaller ds suggests a stronger correlation. 'NC, IC, BC, IC+BC, Whole' mean 
narrow component, intermediate component, broad component, intermediate plus broad component {i.e. narrow component subtracted) and the whole emission line, separately. For [OIII]5007, 'Centre, 
Blue, Whole' mean the central lower velocity shift component, higher velocity shift component (often shows blue shift) and the whole emission line. 'Compton' indicates the soft X-ray Comptonisation 
component. 'Power law' indicates the hard X-ray Comptonisation component. 'Com-|-Pow' indicates the sum of soft and hard X-ray Comptonisation emission (i.e. corona emission). The 'SED Component 
Fraction' indicates the fraction of each component in the bolometric luminosity. '/rac2— iofeev'(i-^- ^2-iOk v) ^^ *^® fraction of 2-10 keV luminosity related to bolometric luminosity. 
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